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ABSTRACT 

We show that if a sample of galaxy clusters is complete above some mass 
threshold, then hierarchical clustering theories for structure formation pre- 
dict its autocorrelation function to be determined purely by the cluster abun- 
dance and by the spectrum of linear density fluctuations. Thus if the shape of 
the initial fluctuation spectrum is known, its amplitude as can be estimated 
directly from the correlation length of a cluster sample in a way which is 
independent of the value of Qo- If the cluster mass corresponding to the 
sample threshold is also known, it provides an independent estimate of the 
quantity (Ts^Iq'^. Thus cluster data should allow both erg and Qq to be deter- 
mined observationally. We explore these questions using N-body simulations 
together with a simple but accurate analytical model based on extensions 
of Press- Schechter theory. Applying our results to currently available data 
we find that if the linear fluctuation spectrum has a shape similar to that 
suggested by the APM galaxy survey, then a correlation length ro in excess 
of 20/i~^Mpc for Abell clusters would require us > 1, while ro < 15/i~^Mpc 
would require as < 0.5. With conventional estimates of the relevant mass 
threshold these imply Qq ^ 0.3 and JIq 1 respectively. 
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1. INTRODUCTION 

Observations have shown that clusters of galaxies are strongly clustered 
in space. The cluster-cluster two-point correlation function, ^cc(^)) is roughly 
a power law, ^cc('") — with a ~ 1.8 and with a correlation length tq 

much larger than that of galaxies (see Bahcall 1988 for a review). The exact 
value of ro is, however, still controversial. Redshift surveys of clusters from 
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the Abell catalogue (Abell 1958) and its southern extension (Abell, Corwin 
& Olowin 1989, AGO) give tq ~ 20-25 h-^Mpc Q (Bahcall & Soneira 1983; 
Klypin & Kopylov 1983; Postman, Huchra & Geller 1992; Peacock & West 
1992). These data also show that the correlation length ro increases with 
cluster richness TZ according to a scaling relation tq ~ OAd, where d is the 
mean intercluster separation and is related to the mean space density n by 
d = (Bahcall & Burgett 1986; Bahcall & West 1992). The rehability 

of these results has, however, been questioned by a number of authors (e.g. 
Sutherland 1988; Dekel et al. 1989; Efstathiou et al. 1992b) who argued that 
they may be affected by projection effects intrinsic to the construction of the 
Abell/ AGO catalogues. These authors claim that the correlation function 
of 7?. > 1 Abell clusters may have been significantly overestimated. Other 
authors (e.g. Bahcall & West 1992; Jing, Plionis & Valdarnini 1992; Peacock 
& West 1992) claim that any bias in ^cc due to projection effects is small. 
Recent estimates of ro based on the APM cluster catalogue (Dalton et al. 
1992, 1994) and the Edinburgh- Durham cluster catalogue (Nichol et al. 1992) 
give ro = 13 — 16 h^^Mpc for clusters with (i ~ 35 — 45 h~^Mpc (compared to 
d ~ 55 /I'^Mpc for TZ > 1 Abell clusters). These results are still consistent 
with the TQ-d relation described above. Since these catalogues are too small 
to give an accurate estimate of ro for richer clusters, the controversy over the 
correlation length of rich Abell clusters is still unresolved. 

Despite the uncertainties in current observational results, the correlation 
function of clusters of galaxies remains one of the most important diagnostics 
for models of structure formation. Indeed, the strong observed correlation of 
clusters on large scales is very difficult to reconcile with the standard cold 
dark matter (GDM) model (White et al. 1987; Dahon et al. 1992; Jing et al. 
^Throughout this paper, we write the Hubble constant as Ho — 100/ikms^^ Mpc^^ 
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1993; Mo, Peacock & Xia 1993; Dalton et al. 1994). In the near future, as new 
data from x-ray observations and large digital sky surveys become available, 
cluster correlation functions will provide much more stringent constraints 
on models. In addition, future observations of cluster x-ray properties and 
of gravitational lensing by clusters will provide improved estimates for the 
masses of clusters. We can then examine ^cc ^ explicit function of cluster 
mass. 

It is clearly important to understand how cosmological models are con- 
strained by such observations. In this context, a simple analytic model for 
^cc and for cluster abundances is particularly desirable. 

The observed high amplitude of ^cc is usually considered to be a result of 
clusters being high peaks in the initial density field (Kaiser 1984). Although 
qualitatively correct, this peak theory does not provide an adequate esti- 
mate for the two-point correlation functions of clusters in N-body simulations 
(Croft & Efstathiou 1994; Mann, Heavens & Peacock 1993). Furthermore, it 
is not clear how to derive a mass function for clusters from the peak theory, 
because peaks of small masses may be contained by those of larger masses. 

An alternative scheme, the Press-Schechter formalism (Press & Schechter 
1974, PS) defines objects (halos) to be virialized structures, identified at the 
present time, which have grown from Gaussian initial density fluctuations. 
Although the original derivation of PS is flawed and recent rederivations are 
still far from rigorous (Bower 1991; Bond et al. 1991), the PS formalism gives 
surprisingly good fits to the mass functions of dark matter halos identified in 
N-body simulations of hierarchical clustering (Efstathiou et al. 1988; White, 
Efstathiou & Prenk 1993; Lacey & Cole 1994). Recently, Mo & White (1995) 
show that the PS formalism and its extensions can also be used to derive a 
model for the spatial correlation of dark matter halos in hierarchical models. 
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They find that the halo-halo two-point correlation functions predicted by this 
model agree surprisingly well with those derived from N-body simulations of 
evolution from scale-free initial power spectra. Such a model can obviously 
help us to understand better the clustering properties of clusters of galaxies. 
A related model has been used by Kashlinsky (1987) to explain the richness 
dependence of cluster correlations. 

In this paper, we present additional tests of the model of Mo & White 
for realistic initial power spectra and in the regime relevant to clusters of 
galaxies. We also discuss how these results can be used to relate observations 
of the large-scale clustering of clusters to the fundamental parameters of 
cosmology. We describe our model in Section 2. The N-body simulations 
used to test our model are described in Section 3. Comparisons of model 
predictions with simulation results are presented in Section 4. In Section 5 
we apply our model to real clusters. A brief discussion of our results is given 
in Section 6. 

2. THE MODEL 

In this paper, we consider models in which the universe is spatially fiat (so 
that fio + = Ij where Qq and Q\ are cosmic density parameters referring 
to matter and to the cosmological constant respectively) and dominated by 
cold dark matter. We also assume that the primordial density field is Gaus- 
sian, with a scale-invariant power spectrum. Neglecting the contribution of 
baryonic matter to Qq, the initial power spectrum of density fluctuations can 
be written as (see Bardeen et al. 1986, equation G3): 

P{k) oc kT\k), (la) 
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with 

T{k) = ^""^V^^"^^^^ fl + 3.89g + (16.1g)2 + {5A6qf + {6.71q)^] (16) 



and 

where, following Efstathiou et al. (1992a), we have introduced a shape pa- 
rameter, r = Qo^) for the power spectrum. 

A useful function equivalent to the power spectrum is the rms mass fluc- 
tuation in a spherical top-hat window of radius R: 

o\R)=Y.P{k)W\R-k), (2) 

k 

where W{R, k) is the Fourier transform of the top-hat window. We normal- 
ize the power spectrum by specifying ag = a{8 h~^Mpc). We define a mass 
parameter Mg as the mean mass contained in a sphere with radius 8 h~^Mpc 
in a universe with the critical density, thus Mg = 5.8 x IO^'^H^^Mq. When 
lengths are written in units of h~^Mpc and masses in units of Mg, all quan- 
tities we are interested in are independent of the Hubble constant. A model 
will then be specified by giving {D,q, F, as)- 

In the PS formalism, dark matter halos are defined as spherically sym- 
metric, virialized clumps of dark matter particles. The mass M of a halo is 
related to its initial comoving radius R (measured in current units) by 

M = -^pR\ (3) 

where p is the mean mass density of the universe at present time. Thus the 

rms mass fluctuation on the scale of the halo is a{R). The fraction of matter 
in halos with mass exceeding M is 



— erfc 



72 



(4) 



6 



where v = 5c/<j{R) and erfc(a;) is the complementary error function. The 
critical linear overdensity for virialization, is chosen to be 5c = 1.69, 
irrespective of Qq, as discussed in White et al. (1993, hereafter WEF). The 
PS formula for the comoving number of halos with mass M is 



^ P d Incr 
n M dM = - - 77-T-r^exp 



S 



2(72 



dM 



(5) 



It follows that the total number density of halos with mass exceeding M is 



, , ^, 3 1 5c d\na 



SI 



2(72 



(27r)3/2 jR a{x) d Inx 

Thus for given a{R) the abundance of halos depends only on R. An Qq- 
dependence enters only through the relation between R and M. 

According to Mo & White (1995), the two-point correlation function of 
dark matter halos with mass M is related to that of the mass, ^m, by 



where 



ar) = [b{M)fUr), (7a) 



6(M) = l + i(i/2-i). (76) 

5c 



This is similar to what one gets from the peak-background split argument 
(Efstathiou et al. 1988; Cole & Kaiser 1989). When halos with a range 
of masses are considered, b{M) in equation (7b) should be replaced by its 
average [weighted by n{M)] over M. Using N-body simulations of scale-free 
spectra. Mo & White show that equation (7) holds even for ^m('^) ,^ 1, as 
long as r is not smaller than the Lagrangian radius R of halos. For rich 
clusters of galaxies, this Lagrangian radius is i? ~ 10/i-^Mpc (WEF). Thus 
we expect equation (7) to be valid for S,cc{r) on r ^ 10/i~^Mpc. On these 
scales, im{r) can be represented sufficiently accurately for our purposes by the 
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Fourier transform of the linear power spectrum. For a given shape of power 
spectrum, the abundance of clusters gives a relation between the threshold 
R and the amplitude cxg (through equation 6). The autocorrelation function 
gives another such relation through the definition p = 6c/(t{R) and equation 
(7). These two relations determine both R and ag. An estimate of the cluster 
mass at threshold then determines 

3. N-BODY SIMULATIONS 

We use four sets of P"^M N-body simulations to test the theories presented 
in the previous section. Each simulation can be characterized by three model 
parameters (r^o^r, erg) (as discussed in Section 2), and three simulation pa- 
rameters: the box size L (in h~^Mpc), the number of simulation particles N 
and the effective force resolution r] (in h~^Mpc) . For the first set of sim- 
ulations, we use (Qo,r,(78) = (0.3,0.225,1) and {L,N,r]) = (400, 100^ 0.2). 
Six realizations were run for this set. The second set consists of five realiza- 
tions of the SCDM universe with (Qq, T, as) = (1.0, 0.5, 1.24) and (L, N, rj) = 
(300, 128^ 0.23). The third is a single realization with ((^o,r) = (1,0.2) 
and (L, N, rj) = (256, 128^, 0.2). Analysis is made for two output times with 
(Tg = 1 and (jg = 0.5. The final set is another single simulation with the 
same initial density spectrum (i.e. with F = 0.2) and the same simulation 
parameters as the third set, but for a low density fiat universe with Qq = 0.2. 
The two output times for this case correspond to erg = 1.07 and 0.5, respec- 
tively. The mass of each individual particle is 5.3 x lO^^M©, 3.6 x IO^^Mq, 
2.2 X and 4.5 x IO^^Mq in these four sets of simulations. 

4. COMPARISON WITH N-BODY SIMULATIONS 

For most of our discussion, we will use clusters identified by the standard 
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friends-of-friends (FOF) group finder. A cluster in a simulation is then de- 
fined as all mass particles that are connected by joining particle pairs with 
separations smaller than a given linkage length, /, in units of the mean inter- 
particle separation. There is no a priori reason for a particular choice of I. 
Following common practice, we choose I — 0.2 (see e.g. Davis et al. 1985; Ef- 
stathiou et al. 1988; Lacey & Cole 1994). The mean density within a cluster 
(S) is then of the order 200, in rough agreement with what one obtains for 
a spherically symmetric, virialized object in an Einstein-de Sitter universe 
(Gunn & Gott 1972). For a spatially fiat universe with f2o = 0.2, S predicted 
by the virialization model is about 2.5 times higher (WEF), and a smaller 
linkage length might seem more reasonable. However, the real virialization 
process must be much more complicated than a spherical accretion model, 
and the appropriate value of / may depend on cluster density profiles, which 
are different in different models (e.g. Jing et al. 1995). In the following 
we will see that our choice of I gives a mass function that is in good agree- 
ment with that predicted by the PS formalism, even for a low-density flat 
universe. For comparison, we will also present some results obtained using 
the spherical-overdensity (SO) grouping algorithm invented by Lacey & Cole. 
This algorithm is based on finding spherical regions with a certain predefined 
mean overdensity k,. A local density near each particle is needed to provide 
an initial list of possible halo centres, and is defined as 3{N + l)/{47rr%), 
where r^ is the distance to the N'th nearest neighbour. Further details may 
be found in Lacey & Cole (1994). We follow them in choosing k — 180 and 
= 10. In this algorithm the mass of a cluster is simply the number of 
particles within the bounding sphere. 
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4.1. Cluster number densities 

In Figure 1 the cumulative number densities of clusters in the simula- 
tions (circles) are compared directly with the predictions of the PS formalism 
(dashed curves). The error bars in Fig. la and lb represent the la standard 
deviations between different realizations, whereas those in Fig.lc and Id rep- 
resent la Poisson fluctuations of the cumulative number density. This figure 
shows that the model prediction agrees with the simulation results reason- 
ably well over a wide range of masses. The abundance range relevant for 
clusters of galaxies is n{> M) = IQ-'^-^-lQ-^ih-^Mpc)-^. The model works 
extremely well for these abundances. This is a nontrivial result, because the 
same 6c and / have been used for all cases. We have tried using / = 0.15 for 
models with Qq = 0.2, and found that the agreement between the resulting 
mass function and model prediction is considerably worse. For comparison, 
we also show in Fig.lc and Id the results for clusters identified by the SO 
group finder (crosses). It is clear that the mass function given by this group 
finder is similar to that given by the FOF group finder with linkage radius 
I — 0.2. This result is in agreement with that obtained by Lacey & Cole 
(1994) for scale- free initial density spectra evolved in an Einstein-de Sitter 
universe. 

4.2. The two-point correlation function 

The two-point correlation function of clusters in the simulations is esti- 
mated directly from pair counts: 

where AP(r) is the average number of neighbors, per cluster, with separation 
in the range r ± Ar/2 and n is the mean number density of clusters in the 
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sample. The circles in Figures 2-7 show ,^cc for the FOF clusters identified 
in the N-body simulations. Results are shown for clusters in different mass 
bins, as indicated in the panels by the number of mass particles, Np. The 
corresponding masses of clusters can be obtained by multiplying by the mass 
of each individual particle. The error bars in Figs. 2 and 3 represent the la 
scatter between realizations, while those in Figs. 4- 7 represents a lo" error 
based on Poisson statistics: A^cc('^) = [AP(r)]^/^/[47mr^Ar]. For compar- 
ison, we show as sohd curves the fits of the data points at r > 10 h~^Mpc 
to a model, ^ccir) = A^ya{r), where A is constant. In this fitting, each data 
point is weighted by the inverse of its bootstrap error, as is discussed in Mo, 
Jing & Borner (1992). The crosses in Figs. 4-7 show ,^cc for SO clusters with 
the same mass as the corresponding FOF clusters shown in the same panel. 
In agreement with the results obtained by Mo & White (1995) for scale-free 
initial density spectra, the correlation functions for FOF and SO clusters of 
the same mass are quite similar. The dashed curves in Figs. 2-7 show the 
predictions of equation (7). The figures show that theory and simulations 
agree remarkably well for all cases. 

5. APPLICATION TO REAL CLUSTERS 

As shown in Section 2, for a given shape of power spectrum the autocor- 
relation function and the number density of clusters with Lagrangian scale 
exceeding R are both determined purely by erg. A comparison of the two 
allows R to be eliminated, determining erg in a way which is independent of 

Figure 8 shows the relation between the correlation length ro and the in- 
tercluster separation d for two sets of linear power spectra. One has F — 0.5 
(Fig. 8a), as in the standard CDM model, the other has F = 0.2 (Fig.8b), 
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which is similar to the shape obtained from the angular two-point correla- 
tion function of galaxies in the APM survey (see Efstathiou, Sutherland & 
Maddox, 1990). It is clear that a stringent constraint on Ug can be obtained 
from an accurate measurement of Tq. The squares are observational results 
for different cluster samples as compiled by Bahcall & Cen (1992), while the 
cross shows the result of a recent analysis by Dalton et al. (1994). The data 
point at d = 94 h'^Mpc is based on Abell TZ > 2 clusters, and is very un- 
certain due to the small size of the sample (Peacock & West 1992). As one 
can see from Fig.Sa, the SCDM model with as — 0.5 does not have enough 
large-scale power to match the observational data. Indeed, if ro > 20 h~^Mpc 
for clusters with d ~ 55 h^^Mpc, then SCDM models with any reasonable ag 
can be ruled out. In contrast, models with F = 0.2 are consistent with such 
correlation lengths provided ag ^ 1 (Fig. 8b). As one can see from Fig.Sb, 
in order to obtain an accurate value of as, it is crucial to measure ro accu- 
rately for such rich clusters. If ro for clusters with d ^ 55 h'^Mpc is actu- 
ally significantly overestimated and its true value is smaller than 15 h~^Mpc 
(Sutherland 1988; Efstathiou et al. 1992b), then as < 0.5 is required for the 
power spectrum shape inferred from the APM galaxy sample. At present, the 
observational data are too uncertain to provide reliable constraints. In the 
future, as new data from digital sky surveys and x-ray observations become 
available, these arguments should provide a determination of as which can 
be compared directly with the amplitude on larger scales inferred from the 
COBE measurements of cosmic microwave background anisotropy (Smoot et 
al. 1992; Wright et al. 1994). 

Based on the observed masses and abundances of rich clusters of galaxies, 
WEF obtained the constraint (Jgilo ~ 0.57, for a spatially flat universe. A 
similar constraint can be obtained from the observed masses and correlation 
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functions. As discussed in WEF, the Lagrangian radius of rich clusters, R, 
is about rg = 8 /i^^Mpc. For R near rg the rms mass fluctuation, a{R), can 
be approximated by 



where the index 7 measures the local slope of the fluctuation spectrum, and is 
given by 7 = 0.68+0.4r for the CDM-like spectra discussed here (Efstathiou, 
Bond & White, 1992a). Using equation (9) and equation (7) we obtain 



where is the mass two-point correlation function for a linear power spec- 
trum with (Tg = 1; A'/g = 5.8 x lO^'^h^^MQ, as introduced in Section 2. For 
rich clusters, the first term in the r.h.s. of equation (10) can be neglected, 
and it follows that 



■M- 


27/3 




[ms\ 







For rich clusters, M ^ Mg (WEF) and (^cc/^n) ~ 10. Wc get agOg-^ ^ 0.53 
for r = 0.2. This result is very similar to that obtained by WEF from cluster 
abundances. 

Combining this result (or the result of WEF) with our result for erg, one 
can in principle get an estimate for Qq- Indeed, if the power spectrum of 
mass density fluctuations has the shape suggested by the APM survey (i.e. 
r 0.2), then ro ^ 20/i~^Mpc for Abell clusters would require JIq ~ 0.3, 
while ro ^ 15h~^Mpc would require Qq ~ 1- With better data on cluster 
masses and correlations, these arguments should provide an estimate of flo- 



a{R) ^ as{rs/Ry, 



(9) 




(10) 
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6. DISCUSSION 



Our simulations show that the mass function and the autocorrelation 
function of clusters can be quite accurately predicted by our Press- Schechter 
model. As a result this model provides a simple and intuitive way to un- 
derstand how cluster data constrain ag and Qq. A significant uncertainty 
remains, however, in the operational definition of a galaxy cluster. In our 
model, as in most theoretical studies, clusters are selected according to mass 
and overdensity thresholds. The two group-finding algorithms we have tested 
give similar mass functions and two-point correlation functions. Unfortu- 
nately, the selection criteria in real observational samples can be far more 
complex than assumed by these algorithms, and may differ from catalogue to 
catalogue. In a recent paper. Eke et al. (1995) used N-body simulations of 
the SCDM cosmogony to examine how the two-point correlation function of 
a cluster sample depends on the way it is selected. They found small but sig- 
nificant variations between cluster samples defined in different ways. Their 
results suggest that selection effects must be considered carefully, when a rig- 
orous comparison between models and observations is made. In the future, 
when cluster samples with well defined selection criteria are available, one 
should use simulations that take into account these selection criteria in as 
much detail as possible in order to calibrate the kind of measurements we are 
suggesting. The observed masses and correlations of rich clusters should then 
provide direct and independent measurements of the density of the Universe 
and of the amplitude of mass fiuctuations. 
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Figure 1. Cumulative mass function of FOF clusters in N-body simulations 
(circles), compared to the predictions of the PS formalism (dashed curves). 
Error bars in (a) and (b) denote the scatter between different realizations; 
those in (c) and (d) denote Poisson errors in the number density. In (c) and 
(d), higher curves correspond to higher values of ag. For comparison, results 
for SO clusters are also ploted as crosses in (c) and (d). 




Figure 2. Two-point correlation functions of clusters of galaxies identified by 
the FOF group finder in N-body simulations (circles), compared to model pre- 
dictions (dashed curves). This is for a model with (Oq) T, ug) = (0.3, 0.225, 1). 
Solid curves show fits of the data points at r > 10 h~^Mpc to a model in which 
^cc(^) is proportional to the two-point correlation function corresponding to 
the linear power spectrum. The error bars represent the la scatter between 
different realizations. The value of rUp denotes the mass of a single mass 
particle in the simulation. Results are shown for clusters in different mass 
ranges, as indicated by the number of particles, Np, contained in these clus- 
ters. The value of d denotes the mean intercluster separation for clusters 
with masses exceeding the lower limit of each mass bin. 




Figure 3. The same as Figure 2, for a model with {Qq, F, erg) = (1, 0.5, 1.24). 




Figure 4. The same as Figure 2, for a model with (fio^Tjas) = (1,0.2, 1). 
Unhke in Fig. 2, error bars here denote Poisson fluctuations. Crosses show 
the same results for SO clusters. 




Figure 5. The same as Figure 4, for a model with [flQ, F, Ug) = (1, 0.2, 0.5). 
Crosses show the same results for SO clusters. 




Figure 6. The same as Figure 4, for a model with (Qq, T, as) = 
(0.2,0.2, 1.07). Crosses show the same results for SO clusters. 




Figure 7. The same as Figure 4, for a model with [flQ, F, Ug) — (0.2, 0.2, 0.5). 
Crosses show the same results for SO clusters. 




Figure 8. Correlation lengtli function of the mean intercluster sep- 

aration d, for two sets of models with T = 0.5 (a) and T = 0.2 (b). For a 
given r, the d-ro relation does not depend on D,q, but only on erg. The data 
points (squares) are observational results compiled by Bahcall & Cen (1992), 
and the error bars are la. In addition, the result of a recent determination 
by Dalton et al. (1994) for APM clusters is also plotted (the cross) together 
with its 2(7 error bar. The dotted lines show the relation tq = OAd. 



